Efficient Quantum Circuits for Schur and Clebsch-Gordan Transforms 
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The Schur basis on n d-dimensional quantum systems is a generahzation of the total angular 
momentum basis that is useful for exploiting symmetry under permutations or collective unitary 
rotations. We present efficient (size poly(n, d, log(l/e)) for accuracy e) quantum circuits for the 
Schur transform, which is the change of basis between the computational and the Schur bases. These 
circuits are based on efficient circuits for the Clebsch-Gordan transformation. We also present an 
efficient circuit for a limited version of the Schur transform in which one needs only to project onto 
different Schur subspaces. This second circuit is based on a generalization of phase estimation to 
any nonabelian finite group for which there exists a fast quantum Fourier transform. 

PACS numbers: 03.67.-a,03.67.Lx,03.67.Mn 



A key component of quantum algorithms is their abil- 
ity to reveal information stored in non-local degrees 
of freedom. In particular, one of the most important 
building blocks known is the quantum Fourier transform 
(QFT)[1], an efficient circuit construction for conversion 
between discrete position and momentum bases. The 
QFT converts a vector of 2" amplitudes in 0{n^) steps, 
in contrast to the 0{n2^) which would be required clas- 
sically. 

Another elementary basis change important in quan- 
tum physics is between independent local states and 
those of definite total angular momentum. When two 
identical spins interact with a global excitation, due to 
their permutation symmetry they appear as a singlet or a 
triplet to the external interaction. Such states of definite 
permutation symmetry can naturally hold entanglement, 
a physical resource central to quantum information. 

The basis transformation defined by permutation sym- 
metry is also central to a plethora of quantum in- 
formation protocols. These include methods to esti- 
mate the spectrum of a density operator [2], achieve op- 
timal quantum hypothesis testing[3], perform univer- 
sal quantum source coding [4], concentrate entanglement 
in a distortion- free manner [5], create decoherence-free 
states [6], and communicate without a shared reference 
frame[7]. This deep connection is a natural conse- 
quence of the underlying states and random variables 
being independent and identically distributed. However, 
unlike the QFT, the complexity of this transform has 
been unknown, rendering protocols which use it non- 
constructive. Also, to be useful in quantum algorithms, 
this transform must be efficient, that is, constructible 
for an n spin system using a quantum circuit of depth 
polynomial (versus exponential) in n. 

Here, we resolve this problem by giving an efficient 
construction for performing a transformation between lo- 
cal and total angular momentum descriptions of a set of 
n d-dimensional systems (n "qudits"), for arbitrary n 



and d. This is achieved using a quantum circuit of size 
poly(n, d, log(l/e)) for accuracy e. We believe that this 
basis change, which we call the Schur transform, is im- 
portant not only for quantum information, but also as a 
new building block for future quantum algorithms. This 
is demonstrated by a close connection between the Schur 
transform and a generalization of the quantum phase esti- 
mation algorithm[8] for non-abelian groups, as discussed 
below. 

The Schur Transform — Consider a system of n qu- 
dits, each with a standard local ("computational") basis 
\i), i — I . . . d. The Schur transform relates transforms 
on the system performed by local d-dimensional unitary 
operations to those performed by permutation of the qu- 
dits. Recall that the symmetric group iS„ is the group 
of all permutations of n objects. This group is naturally 
represented in our system by 

P(7r)|iii2 • --in) = |V-i(l)«7r-i(2) " ' ■«7r-i(n)) , (1) 

where tt G iS„ is a permutation and \i1i2 ■ ■ ■) is shorthand 
for (g) 1*2) (8) . . .. Let Ud denote the group of d x d 
unitary operators. This group is naturally represented in 
our system by 



Q{U)\tii2---ln)^U\ti)(g,U\l2)- 



(S)U\in), (2) 



where U EUd- 

The Schur transform is based on Schur duality, a well 
known [9] and powerful way to relate the representation 
theory of P(7r) and Q(C/). For example, consider the 
case of two qubits [n = 2, d = 2). The two-qubit 
Hilbert space (C^)®^ decomposes under Q into a one- 
dimensional spin-0 singlet space and a three-dimensional 
spin-1 triplet space. Both of these are irreducible rep- 
resentations (irreps) of U2, but they also happen to be 
irreps of 82- The singlet state changes sign under per- 
mutation of the two spins, and the triplet states are in- 
variant under permutation. These correspond to the sign 



2 



Psign and the trivial "Ptriviai irrcps of 52 , and thus we can 
write (C2)®2 ^ (Q^ trivial) © (Qo O Psign), whcrc Qj 
is the spin- J irrep of U2- 

This nmisual coincidence between the two representa- 
tions exists for an arbitrary number of qudits, becom- 
ing quite non-trivial for larger n and d. For example, 
the Hilbert space of three qubits (n = 3, d = 2) de- 
composes into (Q3/2 ® T'triviai) © (Qi/2 ® ^2,i), whcrc 
7^2, 1 denotes a particular two-dimensional mixed symme- 
try irrep of S^. In terms of the original (local) basis 
the Qi/2 (8iP2,i space contains two spin-1/2 objects, one 
spanned by |110) + a.'|011) -I- w*|101) (suppressing nor- 
malization) and |001) -h w|100) w*|010), and the other 
obtained by replacing u) = e^'^^l'^ with w*. These two 
spaces correspond to the two mixed symmetries in 7^2,1 • 

The general theorem of Schur duality states that for 
any (integer) d and n, 
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AePart[n,d] 



(3) 



where A is chosen from the set of possible partitions of n 
into < d parts. Thus, there exists a basis for (C^)®" with 
states |A, (7A,PA)sch where A labels the subspaces Q\®V\ 
and q\ £ Q\ and -px £ V\ label bases for Qx and Vx re- 
spectively. We may represent these states of generalized 
definite total angular momentum and permutation mul- 
tiplicity by vectors |A,<7,p) in the computational basis, 
with bit strings A, q, and p. Note that dim(QA) and 
dim(7'A) vary with A, and so and |p) are padded; this 
requires only constant spatial overhead. 

Just as in the examples above, the Schur basis states 
A,qA,PA)sch are superpositions of the n qudit computa- 
tional basis states . . .in), 



Wq\,P\)sch= ^ [Usch 



|«1«2 • ■ •■*n) ■ (4) 



By the isomorphism of Eq.(3), this defines a unitary 
transformation Usch (with matrix elements as given), the 
Schur transform we desire. Usch maps the computational 
basis to the \X,q,p) representation of the Schur basis. 

Applying the Schur transform extracts A, q, and p 
values for a given state, allowing the values be manip- 
ulated like any other quantum data. For example, since 
maximally entangled states are invariant under permu- 
tation, if \^P)ab is a bipartite partially entangled state, 
then a universal scheme for entanglement concentration 
from 1^)®" is given by both parties performing the Schur 
transform Usch, measuring |A), discarding \q) and noting 
that \p) contains a maximally entangled state of Schmidt 
rank dim('pA)[5]. It is then straightforward to map \p) 
reversibly onto the integers {1, . . . , dim(P;^)} to obtain 
the state Ylf^-^ ^ = dimCPA) [10]. 

The defining property of Usch is that it reduces the 
action of Q and P into irreps. For any n ^ Sn and any 



U e Ud, P(7r) and Q(t/) commute, so we can express 

both reductions at once as 

UschQ(f/)P(7r)U^,h = E l^)(^l®qA(C/)®PA(7r), 

AePart(c(,n) 

(5) 

where qa and pa are irreps of Ud and iS„ respectively. 

Quantum Circuit for the Schur Transform — We con- 
struct a quantum circuit [11] for Usch in two stages, first 
for d = 2, then generalizing to d > 2. Each of these 
constructions follows an iterative structure, in which the 
Schur transform on n qudits is realized using 0(n) ele- 
mentary steps, each of which adds a single qudit to an 
existing Schur state of the form |A, q,p). 

For = 2, this elementary step is familiar from ba- 
sic quantum mechanics, because it involves simple addi- 
tion of angular momentum, following the prescription for 
calculation of the Clebsch-Gordan (CG) cocflicicnts[12]. 
In this case, A and q can be conveniently denoted by 
half-integers J and m (with |m| < J < n/2) which 
give the total angular momentTim and the 2;-component 
of angular momentum respectively. And in terms of 
J, the CG transform takes as input |J, m) and a single 
spin |,s), and outputs a linear combination of the states 
\J ± 1/2, m ± 1/2). The amplitudes of the linear combi- 
nation are readily computed using the usual ladder op- 
erators for raising and lowering angular momenta. In 
addition, however, we must distinguish between multiple 
distinct pathways which add up to give the same total J, 
as demonstrated by the three qubit example above. In 
fact, it is the permutation symmetry of these pathways 
which give rise to Vj, and thus we track the pathway 
with another output label p = J' — J. 
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FIG. 1: Quantum circuit implementing Ucg to convert be- 
tween the I J, m)\s) and \J' ,m' ,p) bases, for the d = 2 (qubit) 
case. Following standard conventions[ll], time goes from left 
to right, the | J) and |m) wires hold multiple qubits, and |s) is 
one qubit. The controlled X operation Cx adds the control 
to the target qubits, i.e. Cx\s)\rn) = \s)\m + s). The doubly 
controlled Ry{0j^m') gate implements the rotation given by 
Eq. (6) using the J and m' qubits. 

Putting this together, we can define an elementary 
Clebsch-Gordan transform step Ucg as a rotation be- 
tween two specific basis states, 

\J'_,m',p=-^) 



\J+,m',p = +\ 

cos 9j^„i' - 
sin Oj^m' 



-sm 6j.,n 
cos Oj^rn' 



J,m+)|s 
lJ,m_)|s 



21 



(6) 
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where 



J ± 1/2, m± = m' ± 1/2, and cosOj^m' = 

UcG can be realized with three gates in a 



I J+m' + 1/2 
V 2,7+1 • 

quantum circuit, as shown in Fig. 1, using as one gate 
a controlled rotation about y by angle 6j^rn'- This angle 
is computed using usual quantum and reversible circuit 
techniques[ll] with error e, using poly(log(l/e)) standard 
circuit elements. 

The full Schur transform is implemented by cascading 
UcG as shown in Fig. 2. The complexity of this circuit 
is thus 0{n ■ polylog(l/e)). In the circuit, p\P2 ■ ■ -Pn is 
an encoding in Young's orthogonal basis [13] for p, the 
label capturing the permutation symmetry of the state. 
This basis is an example of a subgroup adapted basis, a 
useful type of basis used by Beals in the fast quantum 
Fourier transform over the symmetric group [13]. For d = 
2, the total angular momentum of the state, J, gives the 
partition A. 
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FIG. 2: Quantum circuit for the Schur transformation Usch, 
transforming between | ) and I J, m,p). 

Construction of the Schur transform for d > 2 fol- 
lows the same ideas as for = 2, but is complicated by 
the necessary notation; the principle challenge is show- 
ing that the elementary Ucg steps for d > 2 can be 
computed in poly(rf) steps. Usch is constructed clS cl CclS" 
cade of 0(n) Ucg transforms, just as for d = 2. Each 
Ucg combines an arbitrary irrep oi U^, a multi-qudit 
state |A, g), with a single qudit state \ik), to obtain a 
multi-qudit superposition of new irrcps oiUd, |A', q'). Si- 
multaneously, the permutation labels \p) are constructed; 
equivalently we could save the values of A that we gen- 
erate in each step. Ucg can be computed efficiently be- 
cause of a recursive relationship between Uqg for x 
and that of Ud-i x U^-i in terms of reduced Wigner 
coefEcicnts[14]. Crucially, there is an efficient classical 
algorithm for the computation of the reduced Wigner 
coefficients [15] needed for Ucg- Specific details of this 
calculation are given elsewhere [10]. The complexity of 
the full Schur transform is thus found to be polynomial 
in n, d, and log(e~^). 

Circuit for Schur Basis Measurement — Unitary trans- 
forms enable measurements in alternate bases; for such 
applications, simpler quantum circuits can be employed, 
as we now demonstrate for measuring A in the Schur ba- 
sis. Fascinatingly, the quantum circuit we construct uti- 



lizes the same structure as that of the quantum factor- 
ing algorithm[l], and its generalizations to non-abelian 
groups[13, 16-18]. 

Measurement of A projects a state onto the different 
spaces Q\ ®V\, as defined in Eq. (3). For any repre- 
sentation of a group for which there exists a quantum 
Fourier transform, we give a circuit for performing this 
projective measurement onto irreps in Fig. 3; when spe- 
cialized to the symmetric group and the representation 
P(7r), this can be used to measure A. 
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FIG. 3: Quantum circuit used in measurement of an irrep of 
a group G- 

The circuit uses the following gate elements. Let Q 
be an arbitrary finite group over which there exists an 
efficient quantum circuit for the quantum Fourier trans- 
form [16], Fg. Define ^ to be a representative set of 
inequivalent irreps of Q and let dp denote the dimension 
of the irrep p € G. Fg is then a unitary transform from 
a space spanned by group elements, \g), g ^ Q to a space 
spanned by irrep labels and row and column indices, \p, 
a, (3), p & Q, and a,/3 = 1 . . .dp. Specifically, 



Vp{g) \p,a,P){a\ (7) 

-I ap 



where Vp{g) is the entry of the ath row and /Jth 

-I a/3 

column of the irrep p evaluated at e ^. Finally, let 

R-(.9)j ,9 £ be a generic representation of this group Q 
and suppose we can efficiently enact the controlled R((?) 
operation, Cr = Y^g^g I^Xfl ® R'(S')- Fg, its inverse, 
and Cr arc the circuit elements used to construct the 
circuit in Fig. 3. 

This circuit is employed in the following manner, to 
measure an irrep of Q. Consider the action of the circuit 
in Fig. 3 when a generic input state \tj}) is fed into the 
space upon which the generic representation of the group 
R(.g) acts. Since R(,g) is a generic representation of the 
group Q, the space upon which this representation acts 
can be decomposed into diff'erent irreps of Q. Let the 
multiplicity of the /ith irrep (p. & Q) in R(.9) be n^. It is 
possible that an irrep n does not appear at all, in which 
case n„ = 0. Then 



(8) 



= 1 k=l 
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meaning wc can expand the generic input state over 
a basis labeled by the irrep labels ji, multiplicity labels 
for the irrep j, and a label for the spaces upon which the 
irrep acts k. Since this basis fully reduces R(5), 



R-(5)IV') = ^^^c^^,3,k\|J'^3)^^^^{9)\k)K■ 
The output of the circuit in Fig. 3 is thus 



(9) 



\0\ 



vp{g) |p,a,/3) ®R(5)|V). 



(10) 

Using the reducible action of R((?) on |^/;) given by Eq.(9) 
along with the orthogonality relationships for irreps[9], 
Eq.(lO) can be reexpressed as 



o = EEE ■ 

^£0 j=l a,0=l 



|/x,a,/3} (g) |/x,j,/3)R. (11) 



The output \(j)ont) has multiple interesting properties 

which we can now exploit. Measuring the first register 
(the irrep label index) produces outcome ii with probabil- 
ity X^j=i X^fe=i \cii,j,k\^- This is exactly the probability 
we would obtain if we were to measure the irrep label 
of l'^). Remarkably, this is achieved independent of the 
basis in which Cr is implemented. For the Schur basis, 
taking Q = Sn and letting R = P denote the natural 
representation of the symmetric group given by Eq.(l) 
so that fj, = \, the measurement circuit gives exactly 
the probabilities wc would obtain if we were to compute 
UschlV') then measure A. This works for the Schur trans- 
form over not just qubits, but also over qudits. Fur- 
thermore, for large qudit dimension rf, the new circuit is 
exponentially faster than the full Schur transform Usch 
followed by measurement; in terms of d, Usch has size 
poly((i), while this measurement circuit is polylog((i). 

The circuit in Fig. 3 is a generalization of the phase 
estimation circuit introduced by Kitaev[19], which uses a 
QFT over an abelian group (for phase estimation, a,U{l) 
phase is approximated by the cyclic group Cat). In con- 
trast, our circuit is applied to non-abelian groups; and 
in addition to allowing measurement of group irreps, this 
circuit also allows operations to be performed in the irre- 
ducible basis of R(5f), for arbitrary groups Q. This cannot 
be accomplished directly with Usch, which is specialized 
to the permutation and unitary groups. 

Conclusion — ^We have shown how to efficiently perform 
the Schur transform and shown how a generalization of 
Kitaev's circuit can be used to perform a useful group 
representation transformation. Whenever we have iden- 
tical copies of a quantum system or are averaging over 
the diagonal action of the unitary group this first trans- 
form is often essential. Whenever the symmetry of a 
quantum state is described by a finite group, the gen- 
eralized Kitaev transform comes into play. One of the 



most fundamental problems in quantum computation is 
the search for new quantum algorithms. In this respect, 
there are few unitary transforms which have both an ef- 
ficient quantum circuit and interpretations which might 
allow these transforms to be useful in an algorithm. We 
are hopeful that our circuits will be useful exactly be- 
cause they have such clear group representation theory 
interpretations. 



Acknowledgments — This work was partially funded by 
the NSF Center for Bits and Atoms contract CCR- 
0122419 and the NSF Institute for Quantum Information 
under grant number EIA-0086048; AWH was supported 
by the NSA and ARDA under ARO contract DAAD19- 
01-1-06. We thank Nolan Wallach for useful discussions. 



[lo; 
[11 



[12; 
[13: 



[14 

[is: 
[16: 



[17] 

[18 
[19 



P. W. Shor, in Proceedings of the 37th Symposium on the 
Foundations of Computer Science (IEEE, Los Alamitos, 
CA, 1996), pp. 56-65. 

M. Keyl and R. F. Werner, Phys. Rev. A 64, 052311 

(2001) . 

M. Hayashi, J. Phys. A 35, 10759 (2002). 

M. Hayashi and K. Matsumoto, Phys. Rev. A 66, 022311 

(2002) ; M. Hayashi and K. Matsumoto, Quantum Inform. 
Compu. 2, 519 (2002). 

M. Hayashi and K. Matsumoto (2002), quant- 
ph/0209030. 

E. Knill, R. Laflammc, and L. Viola, Phys. Rev. Lett. 
84, 2525 (2000); J. Kcmpc, D. Bacon, D. A. Lidar, and 
K. B. Whaley, Phys. Rev. A 63, 042307 (2001). 
S. Bartlett, T. Rudolph, and R. Spekkens, Phys. Rev. 
Lett. 91, 027901 (2003). 

R. Clevc, A. Ekcrt, C. Macchiavcllo, and M. Mosca, Proc. 
Roy Soc. London Ser. A 454, 339 (1998). 
R. Goodman and N. Wallach, Representations and In- 
variants of the Classical Croups (Cambridge University 
Press, 1998); J. Chen, J. Ping, and F. Want, Group Rep- 
resentation Theory for Physicists (World Scientific, New 
Jersey, 2002). 

D. Bacon, 1. Chuang, and A. Harrow, (2005). 
M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quatum Information (Cambridge University Press, 
New York, 2000). 

J. J. Sakurai, Modem Quantum Mechanics (Addison- 
Wesley, New York, NY, 1994), chap. 3, p. 214. 
R. Beals, in Proceedings of the 29th Annual ACM Sympo- 
sium on Theory of Computing (ACM Press, New York, 
NY, 1997), pp. 48-53. 
J. D. Louck, Amer. J. Phys. 38, 3 (1970). 
L. C. Biedenharn and J. D. Louck, Commun. Math. Phys. 
8, 89 (1968). 

C. Moore, D. Rockmore, and A. Russell , in Proceedings 
of the 15th Annual ACM-SIAM Symposium on Discrete 
Algorithms (SIAM, Philidelphia, PA, 2004), pp.778- 787. 
A. Kitaev, Russ. Math. Surv. 52, 1191 (1995). 
A. Klappenecker and M. Roetteler, quant-ph/0309121. 
A. Kitaev, quant-ph/9511026. 



